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INTRODUCTION 


Water  migration  rates  in  resin  layers,  sufficiently  thin  to  be 

representative  of  the  resin  between  adherends  in  an  adhesive  joint,  were 

studied  with  support  from  the  US  Army  European  Research  Office  under  grant 

number  DA-F.RO-76-G-068.  Both  polyester  and  epoxy  resin  samples  containing 

entrapped  air  bubbles  were  cast  between  glass  microscope  slides  and 

immersed  in  water.  The  times  at  which  droplets  of  water  appeared  inside 

the  air  bubbles,  distributed  at  varying  distances  from  the  edge  of  the 

resin  layer,  were  noted  allowing  the.  rates  of  water  uptake  to  be 

determined.  The  results,  reported  in  a  seminar  presented  by  K  H  G  Ashbee 

at  AMMRC,  Watertown,  Mass  on  23  Julv  1976,  give  diffusion  coefficients  of 
_  a  -I  _  _7  ■>  _1 

!0  cm~s  for  polyesters  and  10  cm~s  for  epoxies;  both  faster  than 

the  published  diffusion  coefficients  in  the  bulk  resin,  by  about  an  order 

of  magnitude. 


During  the  first  year  of  tenure  of  AFOSR-77-3448,  careful  cons idcrat ion 
was  given  to  the  mechanisms  which  might  explain  the  apparent  enhancement 
of  water  migration.  Four  possibilities  are: 

(i)  increased  water  flow  through  channels  or  crazes  associated  with  the 
infrastructure  of  the  interfaeial  region, 

(ii)  stress  enhanced  migration  resulting  from  the  stresses  induced  by 
shrinkage  that  accompanies  curing, 

(iii)  stress  enhanced  migration  caused  by  differential  contraction  on 
cooling  from  the  curing  temperature, 

(iv)  stress  enhanced  migration  effected  by  the  stress  system  produced  by 
resin  swelling  on  water  uptake. 

Evidence  in  support  of  the  first  of  these  hypotheses  has  not  been  found, 
and  recent  studies  involving  microscopic  examination  certainly  have 
revealed  some  characteristic  of  channels  had  they  been  present. 
Similarly,  the  second  possibility  is  a  most  unlikely  mechanism,  for  curing 
takes  place  above  Tg  where  dimensional  changes  will  almost  certainly  be 
taken  up  by  viscoelastic  flow.  Some  stress  system  will  remain  from  the 
third  mechanism  although  partial  relief  of  the  forces  is  anticipated  by 
self  adjustment  of  the  glue  line  thickness,  but  a  more  important  source  of 
stress  is  that  of  the  fourth  mechanism  above.  This  swelling  of  the  resin 
due  to  water  uptake  has  been  shown  to  be  strongly  inhomogeneous,  and  is 
believed  to  be  the  principal  cause  of  stress  fields  which  account  for  the 
significant  changes  in  the  diffusion  rates  in  fibre  reinforced  resins. 
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The  Newton's  rings  experiment  which  has  been  develc pod  to  investigate 
water  diffusion  is  the  most  precise  way  of  investigating  water  migration 
in  stressed  epoxies,  and  a  theoretical  model  has  been  devised  to  probe  the 
physical  mechanisms.  The  classical  diffusion  equation  has  been  modified 
by  introducing  a  stress  induced  flux  term,  and  the  influence  of  internal 
stresses  on  the  activation  energy  is  considered.  The  epoxy  resin,  at 
operating  temperatures  substantially  below  the  glass  transition,  is 
treated  as  purely  elastic. 

The  urgent  need  for  a  fracture  mechanics  test,  which  might  be  used  to 
observe  the  effects  of  water  uptake  on  in  adhesive  joints,  was  made 
evident  at  the  workshop  on  structural  adhesives,  organised  by  AFOSR  and 
held  at  Vought  Corporation,  Grand  Prairie,  Texas  in  December  1977.  Some 
useful  ideas  have  been  developed  and  after  discarding  the  original  test- 
piece,  the  validity  of  the  latest  self-stressed  fracture  mechanics  test- 
piece  has  been  shown.  There  are  still  some  interesting  features  of  the 
stress  system  to  be  analysed,  as  there  are  differences  between  the  theory 
and  practical  considerations.  For  example,  the  theory  requires  a  set  of 
boundary  conditions  which  cannot  be  achieved  in  practice,  but  can  be 
approximated  to  by  the  application  of  a  quite  different  stress  system. 
When  the  test  is  carried  out  the  stress  system  should  revert  to  that  which 
is  used  in  the  theoretical  approach,  but  a  doubt  has  arisen  that  this 
transition  does  in  fact  occur.  This  work  is  however  at  an  advanced  stage 
and  is  continuing  unsupported  by  outside  funds. 

Much  of  the  work  associated  with  water  uptake  causing  resin  swelling,  and 
with  fracture  mechanics  tests  results  in  defects  of  one  kind  or  another, 
and  a  need  arose  to  generate  direct  evidence  for  the  location  and  size  of 
such  defects.  In  view  of  the  progress  made  in  previous  photoelastic 
imaging  experiments  relating  to  a  project  on  the  physical  mechanisms 
responsible  for  the  weathering  of  epoxy  resins  and  CFR  epoxy  resins 
(US  Army  grant  DA-ERO-76-G-068)  a  natural  development  was  to  use 
photoelastic  imaging  of  interference  caused  by  defects.  The  interference 
patterns  are  generated  in  ultrasonic  wave  trains  and  carry  information 
which  allow  cracks,  voids,  and  holes  to  be  located.  By  acoustically 
coupling  a  specimen  to  a  transparent  visualizing  block,  information  about 
defects  in  opaque  specimens  can  be  found.  Attention  has  been  focused  on  a 
computer  model  which  investigates  the  waves  originating  from  a  transducer 
probe  and  those  waves  which  emanate  from  a  defect  or  fault. 

A  series  of  tests  has  proceeded  in  parallel  with  the  computer  model,  but 
as  attempts  have  been  made  to  detect  smaller  and  smaller  defects, 
increased  frequency  requirements  have  led  to  difficulties.  At  frequencies 
above  10  MHz  there  is  a  significant  increase  in  the  attentuation  of  the 
ultrasound,  and  resin  based  specimens  including  those  for  adhesive  joints 
and  film  composites  dissipate  the  stress  waves  rapidly.  With  some  energy 


always  lost  at  the  Interfaces  of  the  two  phases  of  composite  specimens,  it 
is  anticipated  that  the  development  of  ultrasonic  methods  in  non¬ 
destructive  evaluation  will  be  limited. 


2.  NEKTON'S  RINGS  EXPERIMENT 

2.1  EXPERIMENTATION 

The  method  used  to  investigate  swelling  of  the  adhesive  is  to  measure  the 
displacement  normal  to  an  adhesive  layer.  Accurate  measurements  can  be 
achieved  by  observing  Newton's  rings,  perhaps  more  correctly  described  ns 
Fizeau  interference  fringes,  created  by  having  one  of  the  adherends 
sufficiently  thin  to  allow  flexing  and  in  contact  with  a  reference  flat. 
The  specimen  consists  of  a  thoroughly  cleaned  glass  cover  slip  (the  thin 
flexible  adherend)  bonded  to  a  thoroughly  cleaned  rigid  block  of  glass  or 
metal.  With  monochromatic  light,  a  complex  pattern  of  interference  fringes 
is  formed  from  the  cover  slip  and  optical  flat  reflections,  which 
precisely  define  the  topology  of  the  cover  slip  and  therefore  of  the 
underlying  resin.  Should  the  shape  of  the  cover  slip  change,  as  happens 
when  the  resin  swells  during  water  uptake,  the  pattern  of  interference 
fringes  changes. 

The  experimental  arrangement  is  shown  in  Figure  1.  Light  from  a  mercury 
vapour  discharge  lamp  is  directed  through  an  interference  filter,  a 
collimating  lens,  a  half-silvered  mirror,  an  optical  flat,  and  the  cover 
slip.  Interference  between  incident  and  reflected  beams  occurs  within  the 
variable  thickness  gap  located  between  the  cover  slip  and  optical  flat. 
The  interference  pattern  is  photographed  using  the  light  reflected  into  a 
35  mm  camera  by  the  half-silvered  mirror. 

The  pattern  of  interference  fringes  is  determined  by  the  geometry  of  the 
gap  between  cover  slip  and  optical  flat  and,  to  ensure  that  changes  in  the 
pattern  arise  only  from  distortion  of  the  cover  slip  caused  by  resin 
swelling,  it  is  essential  that  each  run  be  carried  out  without  disturbing 
the  specimen/optical  flat  assembly.  All  of  the  components  are  set  up  on  an 
optical  bench.  Figures  2  and  3  show  general  views  of  the  ambient 
temperature  and  high  temperature  (<100°C)  apparatus  respectively. 

Enhancement  of  fringe  contrast  as  a  consequence  of  multiple  beam 
interference  can  be  obtained  by  fully  silvering  the  upper  surface  of  the 
cover  slip  and  partially  silvering  of  the  lower  surface  of  the  optical 
flat.  A  range  of  reflective  coatings  has  now  been  evaluated  from  which  it 
is  evident  that  a  chromium  layer  is  most  effective,  having  sufficient 
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water  resistance  to  withstand  a  full  length  boiling  water  test.  The  main 
disavantages  of  a  chromium  layer  are  that  it  obscures  direct  observations 
of  the  adhesive  and  is  difficult  to  remove  should  it  need  renewing.  An 
alternative  and  very  successful  method  for  securing  multiple  beam 
interference  and  hence  for  improving  fringe  definition  has  been  realised 
by  using  a  wedge  shaped  optical  flat.  With  a  wedge  angle  of  4°  the  image 
of  the  Newton's  rings  does  not  coincide  with  reflections  from  other  glass 
or  water  surfaces.  The  wedge  is  made  from  flint  glass  with  refractive 
index  1.65  and  has  been  incorporated  into  a  mark  II  version  of  the  hot 
water  apparatus  designed  to  permit  undisturbed  tests  at  any  fixed 
temperature  between  ambient  and  100°C.  This  is  shown  in  Figure  4. 

Whilst  direct  observation  of  this  interference  pattern  can  be  used  to 
measure  increments  of  swelling,  the  changing  geometry  is  obtained  directly 
by  S"perimposing  successive  images  of  the  pattern  photographed  during 
water  uptake,  on  the  pattern  photographed  before  swelling. 
Suporimposition  produces  Moire'  fringes,  the  development  of  which 
faithfully  follows  any  changes  in  the  shape  of  the  cover  slip. 

Moire  fringes  are  generated  when  patterns  with  nearly  identical  periodic 
structures  are  made  to  overlap.  In  this  instance,  photographic  images  of 
the  interference  pattern  form  the  periodic  structures  and  differences  in 
the  pattern,  seen  on  comparing  the  image  photographed  before  swelling  with 
images  taken  during  swelling,  gives  rise  to  the  Moire  effect.  In  Figure 
2,  the  Moire  fringes  are  the  circumferential  lines  running  around  the  rim 
of  the  specimen  and  are  superimposed  upon  the  interference  pattern.  (It 
should  be  noted  that  the  broad  dark  almost  parallel  bands  running,  across 
the  specimen  are  due  to  tilting  with  respect  to  the  optical  flat).  As 
water  diffuses  Into  the  specimen,  new  interference  fringes  are  generated 
at  the  edge  and  migrate  towards  the  centre,  each  fringe  representing  the 
locus  of  points  of  constant  displacement  normal  to  the  joint.  Similarly  in 
the  Moire  pattern,  a  now  Moire  fringe  is  generated  and  migrates  in  from 
the  edge  of  the  specimen  whenever  a  new  interference  fringe  is  generated. 
Thus  successive  Moire  fringes  map  out  the  loci  of  positions  which  differ 
in  displacement  normal  to  the  joint  by  half  a  wavelength.  Hence  the  number 
of  Moire  fringes  N  to  have  passed  a  given  point  is  simply  related  to  the 
normal  displacement  of  the  joint  w  by: 

w  =  N  X 

2  P 


where  A  =  wavelength  of  light 

jj  =  refractive  index  of  the  immersion  liquid 


Specimens  usually  took  the  form  of  a  circular  glass  cover  slip  bonded  to  a 
rigid  block  of  metal  such  as  stainless  steel  or  anodised  titanium. 
However,  some  specimens  were  manufactured  using  cover  slips  which  were 
both  square  in  outline  and  of  different  thickness.  This  was  done  in  order 
to  investigate  the  effects  of  varying  the  mechanical  constraint 
attributable  to  cover  slip  geometry. 


2.2  WATER  ABSORPTION  MODEL 

Consider  the  adhesive  specimen  of  the  previous  Section  in  cylindrical 
polar  coordinates  and  assume  axial  symmetry. 

Let  JJr,t)  denote  the  flux  of  water  molecules  passing  perpendicular 
through  a  reference  surface  of  unit  area  during  unit  time.  may  be 
expressed  in  terms  of  the  local  mass  velocity  of  the  water  molecule  by 

J_(r ,  t )  =  c(r,t).jj(r,t) 

where  c(r,t)  is  the  local  concentration  of  water  molecules,  at  time  t.  No 
mechanism  of  migration  need  be  specified  at  this  point. 

The  aim  is  to  investigate  the  influence  of  the  internal  stress  field  which 
causes  changes  in  both  c(r,t)  and  ujr,t). 

Fick's  second  law, 

5c  =  5  {D  (c)5c  } 

5t  dx  dx 

gives  the  concentration  of  a  diffusing  species  in  one  dimension.  D,  the 
diffusion  coefficient  can  often  be  taken  as  constant,  but  in  the  diffusion 
of  polymers  it  is  a  function  of  concentration.  When  extensive  swelling,  of 
a  polymer  is  caused  by  the  penetrant,  non-Fickian  diffusion  has  been 
observed,  Kwei  and  Zupko  (1969),  Alfrey  et  al  (1966). 

Non-Fickian  behaviour  may  be  directly  related  to  the  influence  of  the 
changing  polymer  structure  on  the  solubility  and  differential  mobility,  or 
may  result  from  the  Internal  stresses  exerted  by  one  part  of  the  medium  on 
another  as  diffusion  proceeds.  Polymers  usually  have  a  wide  spectrum  of 
relaxation  times  associated  with  structural  changes  and  a  sorption  process 
will  be  influenced  by  those  segmental  motions  which  occur  at  about  the 
same  rate  or  slower  than  the  motivating  diffusion  process.  Alfrey,  Gurnee 
and  Lloyd  (1966)  proposed  a  useful  classification  according  to  the 
relative  rates  of  diffusion  and  polymer  relaxation.  Three  classes  are 
distinguished: 
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(I)  Fickian  diffusion  in  which  the  rate  of  diffusion  is  much  less  than 

that  of  relaxation. 

(II)  Diffusion  in  which  the  rate  is  very  rapid  compared  with  the 

relaxation  process. 

(III)  Non-Fickian  or  anomalous  diffusion  which  occurs  when  the  diffusion 
and  relaxation  rates  are  comparable. 

The  following  assumptions  are  made  in  the  development: 

Epoxy  resins  are  assumed  to  be  isotropic  and  homogeneous. 

The  specimen  is  treated  as  having  axial  symmetry  at  all  times. 

Activated  diffusion  forms  the  main  contribution  to  the  migration  of  water 
molecules  within  the  epoxy. 

No  increased  migration  of  water  takes  place  along  the  glass/epoxy 
interface  due  to  any  deterioration  of  interfacial  bonding. 

Water  uptake  produces  stresses  which  are  not  sufficiently  large  to  form 
microcracks  and  channels. 

Effects  of  gravity  are  negligible. 

The  epoxy  resin  behaves  as  perfectly  elastic  for  the  internal  stresses  and 
strains  caused  by  water  uptake. 

The  elastic  constants  are  unchanged  by  the  water  concentrations 
encountered  in  the  epoxy  resin. 

The  compressive  forces  at  the  rim  of  the  specimen  resulting  from, 
inhomogeneous  swelling  have  negligible  effect  on  the  saturation 
concentration  at  the  outer  surface,  and  that  this  equilibrium  saturation 
concentration  is  reached  instantaneously. 

In  the  following,  two  different  models  for  the  effect  of  stress  on  water 
uptake  will  be  considered.  It  will  become  apparent  that  a  very  important 
quantity  in  suet)  models  is  the  concentration  gradient.  A  large 
concentration  gradient  implies  a  large  flux  of  water  molecules  by  Fick's 
First  Law,  which  gives 

J  =  -D  d_c  2.  I 

dx 

In  addition,  however,  a  large  concentration  gradient  implies  a  large 
stress-field  and  thus  possibly  a  significant  stress-induced  flux.  In 
subsequent  paragraphs  consideration  will  be  given  to  the  effect  of  the 
stress-field  on  the  mobility  and  thus  on  the  diffusion  coefficient  D.  A 
stress-induced  flux  term  in  Fick's  First  Law  will  also  be  introduced. 

It  will  be  convenient  to  use  the  word  "system"  to  refer  to  the  epoxy  resin 
contained  between  the  cover-slip  and  the  microscope  slide  plus  any  water 
molecules  tiiat  have  diffused  into  it.  The  word  "specimen"  will  refer  to 
the  complete  model  adhesive  joint,  i.e.  the  epoxy  resin  plus  the  glass 
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cover-slip  and  microscope  slide. 

Now  consider  a  water  molecule  within  the  epoxy  resin.  On  average  it  will 
make  jumps  per  unit  time  with  a  distance  6  per  jump.  Then  the  diffusion 
coefficient  will  be  given  by  an  expression  similar  to 


I)  =  1 
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<t>6 


Now  5  tlie  jump  distance  will  clearly  depend  on  the  state  of  the  system; 
that  is,  one  would  expect  the  jump  distance  in  a  rep i on  under  compression 
to  be  less  than  that  in  an  unstressed  rep, ion.  It  will  be  assumed  for  the 
moment,  however,  that  the  jump  distance  remains  constant  and  cons  iderat i on 
will  bo  piven  to  the  effect  of  stress  on  the  number  of  jumps  per  unit 
time,  <J> . 


Lot  v^  denote  that  volume  within  the  system  which  is  effectively  free  to  a 
diffusinp  species  to  occupy  without  the  necessity  to  perform  work  on  the 
surroundinp  polymer  molecules.  For  an  unstressed  specimen  one  would  expect 
this  volume  to  vary  throuphout  the  system  as  a  result  of  the  fluctuations 
in  thermal  enerpy  of  the  constituent  polymer  molecules.  Therefore , 
consider  the  mean  value  of  v^- ,  to  be  denoted  by  \y ,  which  will  he  referred 
to  as  the  "free  volume".  In  a  repion  of  the  polymer  which  is  under 
compression  would  be  expected  to  decrease  from  its  unstressed  value  and 
in  a  repion  under  tension  one  would  anticipate  Vj-  to  increase.  Thus  for  a 
system  in  which  there  exists  a  stress-field 


vf  =  vf  (r , z , t ) 

Now  returning  to  <f  ,  the  number  of  jumps  per  unit  time,  one  can  see  that 


❖ 


f<EthV 


where  F.^ 


denotes  the  thermal 


molecule  being 


t|(  uriivu.cn  mi-  Liirn-uix  energy  of  the  water 
considered  and  f  Is  a  function  to  be  determined.  If  the  effective  volume 
of  the  water  molecule,  -v  ,  in  the  system  is  greater  than  the  free  volume 
Vj  then  work  will  have  to  be  done  on  the  polymer  molecules  surrounding  the 
free  volume  to  permit  the  water  molecule  to  jump.  The  enerpy  to  perform 
this  work  must  come  from  the  thermal  fluctuations.  The  probability  that 
the  water  molecule  has  thermal  energy  1-1'  ,  is  proportional  to  the 


tit 


— F '  /KT  '  th  ^ 

lloltzmann  factor  e  '  th'*  •  Now  for  a  jump  to  occur  one  must  have  b 

K*  where  K*  denotes  the  energy  equivalent  of  the  minimum  performed  work  in 

jumping  plus,  possibly,  some  constant  energy  term.  Hence  the  probability 

of  jump  Pr  is  proportional  to 


7 


1/ 


e-Eth/kT  dE- 


CO 


E* 


til 


}/{/ 


e-tth/kT  dE' 


th  1 


giving  Pr  “  e~E  where  E*  is  called  the  activation  energy  of  the  process. 
Hence  <t>  01  e  "  '  and,  subsequently,  from  2.2  one  obtains  D  “  e!*  /K  . 

If  one  assumes  that  variations  in  the  jump  distance  6  have  a  negligible 
effect  on  D  compared  to  variations  in  <i>  then  a  possible  expression  for  the 
diffusion  coefficient  would  be 


D 


Doe 


-E*/kT 


2.3 


An  alternative  way  of  analysing  the  diffusional  jump  is  to  think  of  the 
fluctuations  in  local  free  volume  v^  due  to  the  thermal  fluctuations  of 
the  energy  of  the  polymer  chains.  Then  a  water  molecule  will  only  jump  if 
a  sufficiently  large  free  volume  v^  (=  vy)  exists  adjacent  to  its  present 
site.  Thinking  in  terms  of  the  free  volume  as  consisting  of  discrete, 
independent  entities  one  may  use  the  term  hole  instead  of  free  volume.  If 
an  energy  E^  is  associated  with  a  hole  of  size  v^  then  the  probability 
that  a  hole  of  size  v^  can  be  found  is,  according  to  the  Boltzmann 
distribution,  proportional  to  expf-E^/kT).  The  average  size  of  a  hole  in 
the  material  is  then  given  by  the  following  expression 

vf  •  a;  vfe-Ei.«  dvti/{j[WkT  dvti 


Also  the  fraction  of  holes  having  volume  vw  or  greater  is 

{revkTdvf}/{re~Eh/kTdvf} 

L  O  *• 

In  writing  equation  2.4  one  is  assuming  that  the  holes  are  discrete  and 
independent  but  in  reality  this  is  not  so.  Even  if  an  individual  hole  may 
not  be  large  enough  to  accommodate  a  diffusing  molecule,  the  co-operative 
motion  of  several  neighbouring  molecules  may  allow  two  or  more  holes  to 
merge  into  one  hole  large  enough  for  a  diffusional  jump  to  occur.  Tims  the 
fluctuations  in  free  volume  will  be  brought  about  by  complex 
configurational  changes  of  the  polymer  molecules  and  one  would  anticipate 
that  an  expression  involving  the  various  degrees  of  freedom  available  to 
the  polymer  molecules  would  be  obtained  instead  of  2.3.  This  is  in  fact 
true  and  many  expressions  have  been  proposed  (see,  for  example,  Crank  and 
Park  [1968]  and  Barrer  [1941]). 


2.4 


Returning  to  equation  2.3  it  is  assumed  that  E*  can  be  written  in  the  form 


E*  =  E  +  E 
w  c 


2.5 


where  Ew  is  th c  energy  equivalent  of  the  work  done  by  a  water  molecule  in 
a  diffusional  jump  and  Ec  is  some  constant  energy  term  which  includes  that 
part  of  the  energy  required  by  the  water  molecule  to  change  locations  even 
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if  the  free  volume  is  equal  to  or  greater  than  the  volume  v  .  Consider  a 
spherical  cavity  of  radius  r^  at  the  centre  of  a  sphere  of  large  radius  R. 

Let  r^  and  rw  be  given  by 

vf  =  4  rcrf 3,  vw  =  4  uw3 


It  will  be  assumed  that  a  suitable  approximation  for  F.  is 

Ew  =  Prf  -  rw  (vw  -  Vf)  2-6 

where  Pr^  rw  is  the  pressure  required  to  change  the  radius  of  the 

central  cavity  from  r^  to  rw«  Using  elasticity  theory  one  obtains 

Prf  >  rw  =  2 T _ {n  (r , t)  -  1)  2.7 

( 1+v ) 

where  q  =  r  /rf  >  !• 

Now  clearly  as  water  molecules  migrate  through  tiie  system,  and  a 
continuously  changing  stress-field  is  set  up,  then  the  free  volume  v^-  will 
change  from  its  initial  unstressed  value  v ^ . .  In  order  to  estimate  this 
etiange  in  free  volume  v^  consider  a  thin  spherical  shell  of  thickness  s 
where  v^-  is  represented  by  the  spherical  interior  of  this  shell. 

Assuming  that  r^  is  small  then  l /30j^  (r , z , t )  is  approximately  the  average 
value  of  the  normal  component  of  stress  over  the  surface  of  the  shell  with 
centre  at  a  point  (r,z)  at  time  t.  Tims  it  is  reasonable  to  take  the 
pressure  Pg  at  the  surface  to  be  given  by 


Ps(r,z,t)  =  -l/3oii(r,z,t) 

where  is  assumed  known,  and  determined  by  elasticity  theory.  To  remove 

the  z-dependence  of  this  pressure  consider  an  average  over  the  z- 
coordinate : 

h 

P  (r,t)  =  -l/Sh/o,, (r,z,t)dz 


where  0^  =  orr  +  Ogg  +  ozz  and  h  is  the  thickness  of  the  system.  Assuming, 
tiie  thickness  of  the  shell  to  be  small  relative  to  r^  elasticity  theory 
yields 


u  =  -  Ps  (1  -  v)  r{ 
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Thus  the  new  free  volume  v(  „  is  obtained  from  the  old  free  volume  vr  bv 

In  f  o  J 

vfn  -  4  (  (3vfo>1/3  +  u)3 

3  4n 

=  vfQ  (l  -  X(vfn)1/3p<  ) 

S  2.8 

where  X  =  3(J_)^3  ( 1  -  v ) 

4  n  2E 

Thus  combining  2.3,  2.5,  2.7  and  2.8  one  obtains 

D  =  D'Qexp  [~a(n  -  l){v..  -  vfo  +  Xf.Ps>] 

kT  '  2.9 

when  v  >  vc  ,  where 
w  t  n 

D'o  =  D0cxp[-Ec/kT],  n  =  r^  ,  a  =  _2F_,  T>  “  ^’fo>4/3»  vfo  =  4/3nrfo3, 
h  rw  (1+v)  S 

and  Pg  =  -l/3h/a^  dz 
o 

If  vw  <  v.  then  D  =  DQexp [-Ec /kT ]  and  no  swelling  should  be  observed. 

Having  investigated  the  effect  of  stress  on  the  mobility  and  hence  the 
diffusion  coefficient,  the  possibility  of  a  st ress-i nduced  flux  will  now 
bo  considered.  It  should  be  noted  that  these  two  models  proposed,  are 

different  ways  of  modelling  the  same  effect,  i.e.  the  effect  of  the 

internal  stress-field,  resulting  from  water  uptake,  on  the  subsequent 
water  migration.  An  attempt  to  combine  the  two  models  would  be  incorrect. 

Firstly  irreversible  thermodynamics  will  be  used  to  derive  an  equation  for 
the  migrtion  of  water  molecules  which  incorporates  a  term  which  can  be 
interpreted  as  a  stress-induced  flux.  The  equation  to  be  derived  is 

similar  to  that  proposed  by  Frisch  et  al  (1969)  to  model  Case  II  swelling 

(see  Page  3). 

Consider  an  element  of  volume  V  of  the  system.  Let  there  by  n  molecules  of 
water  in  this  volume,  with  velocity  v^  and  denote  their  centre  of  mass 
velocity  by  Vy. 

Then  v..  °  )  v, 

i=l-i 

Let  the  centre  of  mass  velocity  of  the  polymer  molecules  be  v^ .  Then  the 
centre  of  mass  velocity  vc  of  all  molecules  within  the  volume  V  will  be 
given  by 
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o  V  -  n  V  +  d  V 
^  -c  K  W“W  P  p-p 


where  p,  and  p  are  the  densities  of  the  water  molecules  and  the  polvmer 
h  w  h  p  ' 

molecules  respectively  and  p  =  pw  +  p^.  Hie  flux  of  water  molecules 
defined  with  respect  to  the  centre  of  mass  frame  is  pi von  by 


=  cw(Xw 


?c> 


(  = 


-V 


where  c,  (the  concentration  of  water  molecules) 
w 

conservation  of  mass  gives 


p  /p,  and  the 


pPc ■  +  V.(pJ)..  =0  i  =  w,p 

Dt  2.10 


where  f>  =  VX  +  V  .  V 
Dt  6t 


It  will  be  assumed  that  the  equation  of  motion  may  be  written  in  the 
f  orm:  - 


rLliYc)1 

Dt 


=  ^  + 


/  Pi,  1 

k=w,p 


^J2i  j 


Ci ,  j 


1,2,3) 


by  analogy  with  fluid  dynamics,  where  is  the  stress  tensor,  and  the 
forces  Fk  will  be  assumed  zero.  The  energy  equation  for  u,  the  energy  per 
unit  mass,  with  the  exclusion  of  the  centre  of  mass  kinetic  energy  is 


pH  (i(V2  +  u)  =  iLCPii^ii) 

Dt  2  dxj 

o 

Elimination  of  ( vc )  using  the  equation  of  motion  gives 

PDn  =  P±  i^(V<.)i  (i’j  =  1»2*3)  2.11 

Dt  dx  j 

The  Second  Law  of  Thermodynamics  gives 

T  PS  =  DU  +  pPV  -  T  pR  DMk 

Dt  Dt  Dt  k=w,p  Dt  2.1 


where  (i^  and  M  are  the  chemical  potential  and  total 
molecules,  respectively;  and  are  the  chemical 

mass  of  the  polymer  molecules,  respectively;  S  is  the 
the  total  energy.  The  quantity,  p,  in  2.12,  denoting 
defined  to  be  equal  to  -I/3pj^,  and  for  convenience  let 


mass  of  the  water 
potential  and  total 
total  entropy  and  W 
the  pressure  can  be 
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2.13 


Pij  "  'P&ij  +  dij 

where  d ^ ^  =  0. 

Changing  to  intensive  quantities  2.12  becomes 

T  Ds.  *=  Du  +  pDv  -  Xgk  Dck 
Dt  Dt  Dt  k=w,p  Dt 

Now  substituting  2.13  into  2.11  gives 

pDu  =  -p  V.V^,  +  dH.  d(v„)4  (i,j  =  1,2,3) 

Dt 

But  Dp  = 

Dt 

thus  _Du  = 

Dt 

Introducing  2.10  and  2.15  into  2.14  gives 


dxj 


-v~2  Dv  =  -pv.y^ 
Dt 


-p 


Dv  + 
Dt 


iij 

p 


LL^)i 

Sx 


2.14 


2.15 


pT  Ds  =  dij  +  S  PkV.(pJk)  (i,j  =  1,2,3)  2.16 

Dt  9xj  k=w,p 

Now  diffusion  is  an  irreversible  phenomenon  which  can  he  expressed  to  a 
first  approximation  by  linear  phenomenological  relations  of  the  general 
type 

m 

=  T,  L^k  Xk  (i  =  1,2,  ...,  m) 

k=l 


for  an  m-component  system  where  denotes  the  flux  of  the  i1'1  component, 
denotes  a  "force"  (e.g.  concentration  gradient,  potential  gradient), 
and  the  coefficients  Lik  (i,k  =  1,2,..., m)  are  called  the 

phenomenological  coefficients  (see,  for  example,  Bearman  and  Kirkwood 
(1958),  Bearman  (1958),  De  Groot  (1951). 

Now  2.16  can  be  rewritten  in  the  form 

p  Ds  =  X  V.  (UkPik)  +  dij  9(V<;)1  +  T.  dk-Xk 

Dt  k=w,p  T  T  dxy  k=w,p  T  2.17 

where  =  -pVpk  2.18 

for  isothermal  diffusion. 


12 


Note  the  tern  d ^  ■  ^(v^j)  will  hove  no  influence  on  the  form  of  the 

T  ftxj 

phenomenological  relations 


J 

=  L  X 

+ 

L  X 

— w 

ww— w 

wp^-p 

J 

=  L  x 

+ 

L  X 

-p 

pu— w 

PP-P 

since  the  force  in  this  term  (i.e.  A(V  )./c)x.)  has  a  tensorial  character 

- C  1  J 

and  cannot  give  rise  to  vectorial  fluxes  (Curie's  theorem). 

Now  it  can  be  shown  (See  De  Groot  f 1 95 1 J  )  that 


L 

WW 

wp 

L 

ww 

pw 

0, 

0, 


+  I. 

r 

+  L 


;  pp 
Jpp 


o, 

0. 


so  that 


J 

— w 


-J  = 
~P 


Lww(*w  -  V 


2.  IP 


Alternatively,  since  it  was  already  known  that  d_  =  one  could  have 
used  this  in  deriving  2.17  to  give  =  -f>V  (ji  w  -  p  ),X  =  0  yielding  the 
same  results.  Using  2.18  and  2.19  one  obtains 


iw 


Lww("P^w  + 

LwwP  [(  (Up  -  MW»  VP  +  (  <L(,,p 
9  p  T,C  0(1, 


»*w» 


VC  1 
w 

T,P  2. 


Now  C  du  +  C  du  =  0  at  constant  T  and  p  (  -  the  G i bbs-Puhon  relation) 
w  r  w  P  P 

also  C  +  C  =1.  Now  it  is  expected  that  C  >>  C  .  tlien 
w  p  p  w 


also 


(  L 

dci 


(pp  -  pw)) 


=  -  (i 


tip 


-<V 

Jw 


dC, .  T^p 


)( 


c 

-it 

c  <yc 

p  w 


1* 


(^Pp)  «  («v„) 

0  T.C  Op  T,C 

p  1  w  1  l  w 
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tt.1.  i  A 


Thus  2.20  becomes 


iv  =  "  LwwP[(  3ffw)  +  (^w)  V(U  2.21 

aP  Ticw  0PW  TlP 

t- Vi 

For  non-ideal  systems  the  chemical  potential  of  the  ic  component  takes 
the  form 


% 


Pi  =  C.(p,T)  +  RT  log  f .C.  2.22 

where  C^(p,T)  is  independent  of  composition,  CL  denotes  the  concentration 
and  fi  is  the  activity  coefficient,  due  to  0  N  Lewis  (see  for  example, 
Prigogine  [1967]).  Substituting  2.22  into  2.21  one  obtains 

iv  -  -  LwvP  [(  *£*)•  VP  +  II  VCwl 

ep  T  Cw 

Now  Fick's  First  Law  gives 

J  =  -  DVC 

therefore  let  the  diffusion  coefficient  be  defined  by 

D  =  LwPRT  2.23 

Cw 

Then  the  flux  of  water  molecules  in  the  epoxy  resin  is  given  by 

Iw  =  "D<VCw  +  Ek  Vp) 

RT  5p  T  2.24 

where  the  term  (SC) 

3p  T 

will  depend  on  the  interaction  forces  between  the  water  molecules  and 
their  neighbouring  polymer  molecules. 

An  alternative  method  of  deriving  an  equation  of  the  form  2.24  is  to 
consider  the  probability  of  a  water  molecule  jumping  to  the  right  or  to 
the  left  (in  one-dimension).  For  simplicity  consider  the  one-dimensional 
case  only,  with  position  co-ordinate  x.  Suppose  that  the  pressure  in  the 
epoxy  resin  p(=  -  l/3cs^  =  pg)  decreases  with  increasing  x  -  co-ordinate. 
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Then  the  probability  of  a  water  molecule  jumping  to  the  right,  to  a  site 
at  x  +Ax,  is  given  by  an  expression  of  the  form 


PR  =  Aexp^“^al  +  a2P^x  +Ax)}/kT] 

using  equation  2.9.  Similarly  the  probability  of  a  water  molecule  jumping 
to  the  left,  to  a  site  at  x  -  Ax,  is  given  by  an  expression  of  the  form 

PL  =  AexpMnl  +  a2p(x  "  Ax)}/kT] 

Thus  the  flux  of  water  molecules  due  to  the  stress-field  is  given  by 
JS  =  P<PR  -  PL> 

^  pAn0  {p(x  -  Ax)  -  p(x  +  Ax)> 
kT 

~  -2pAn ,  6  Ap 
kT“  5x 

where  Ax  has  been  replaced  by  6  the  mean  jump  distance  and  p  is  some 
suitable  proportionality  factor. 

Then  the  total  flux,  denoted  by  will  be  of  the  form 

JTOT  =  ~D  <ll£  +  ac 

3x  dx 

I 

similar  to  2.24 
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2.3 


RESULTS  AND  DISCUSSION 


After  the  first  two  years  of  exper incntal  work  the  models  described  in  the 
previous  Section  were  set  up.  These  two  models  of  the  possible  effect  of 
an  internal  stress  field  on  water  migration  in  epoxy  resins,  can  be 
represented  by  the  following  two  equations: 

-{aj+a^pfr.t) }/kT 
be  =  V.[nQe  “  VC] 

3t 

and  b£  -  V.[D(VC  +  acVp)] 

at 

where  aj  and  a0  are  constants  with  a  assumed  constant.  For  both  models  it 
is  important  to  have  a  knowledge  of  the  pressure  distribution  p(r,t)  which 
may  bo  obtained  from  elasticity  theory.  Oj-fr.t)  is  most  convenient lv 
arrived  at  by  applying  numerical  methods  such  as  finite  difference  or 
finite  element  methods.  Thus  one  can  proceed  to  values  of  p(r,t)  ^  -l/3e.. 
from  which  the  values  of  tiie  constants  can  be  determined  by  matching 
computed  concentrations,  for  a  given  time,  to  experimentally  evaluated 
concentrations. 

It  is  possible  as  a  first  approximation  to  use  the  normal  stresses 
necessary  to  return  the  cover  slip  to  its  initial  shape  as  the  pressure, 
and  here  we  could  apply  A  E  11  Love's  equation 

p  =  RV4cj 

where  u>  =  6(0)  is  the  axial  displacement,  p  the  normal  stress  (pressure) 
and  R  the  flexural  rigidity.  The  cover  slip  is  treated  as  a  thin  plate. 

The  exact  nature  of  the  effect  of  the  cover  slip  on  the  stress  field  is 
however  a  major  consideration,  since  replacement  of  the  cover  slip  with 
another  microscope  slide  would  allow  little  deformation  and  surface 
stresses  would  become  more  significant. 

The  effect  of  the  cover  slip  is  not  at  all  obvious  as  it  is  chosen  for  its 
flexibility  but  the  tangential  stresses  exerted  by  the  cover  slip  could 
form  an  important  constraint  on  the  system. 

The  first  long  term  test  was  made  at  room  temperature  on  a  joint  between  a 
cover  slip  and  cleaned,  but  not  anodised.  SIC  aluminium.  A  general 
purpose  epoxy  resin  containing  entrapped  air  bubbles  was  used  as  adhesive. 
Hereafter,  this  adhesive  will  be  referred  to  as  adhesive  A.  Its  full 
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identity,  including  hardener  and  acee'erator,  and  the  curing  schedule  used 
for  it,  are  given  in  Table  1. 


RFSiN’  SYSTF.M 

PROPORTIONS 

CHI. 

CURK 

p  .  b .  w . 

Ciba  Geigv 

MY  750  (Diglycidyl  ether 

of  bisphenol  A) 

100 

2  hours 

4  hours 

A 

HY  917  (Methyl  tetrahydro- 

at 

at 

phthalic  anhydride) 

85 

1 00°c 

1 50°C 

DY  062  (Trinmyl  an’monium 

1 

phenate) 

2 

i 

1 

Bloomingdale 

1  hour 

FM  1000  (Nylon  mcdified 

at  i 

B 

epoxy  film) 

r_> 

G 

O 

2  hours 

I  hour 

BR  1009/49  (Phenolic  primer) 

at 

a  t 

20°C 

70°C 

C 

Ciba  Geigy  (modified  epoxy  film) 

1/2 

hour  at  120°c 

Redux  312/5 

n 

Bloomingdale 

l 

hour  at  1  20°C 

FM  73M  (modified  epoxy  film) 

Table  !•  Adhesives  and  Curing  Schedules. 


Later  tests  were  carried  out  on  joints  between  a  cover  slip  and  substrates 
of  both  stainless  steel  and  anodized  titanium  using  the  commerinlly 
available  adhesive  film  system  designated  B  in  Table  1.  All  metal 
surfaces  were  cleaned  and,  where  applicable,  treated  in  accordance  with 
commercially  accepted  practice.  The  metal  based  specimens  were  supplied 
by  Westland  Helicopters  Ltd  in  the  form  of  cleaned  and  pre-prined  coupon 
specimens . 

For  tests  employing  the  adhesivr  film  systems  designated  C  and  D  in  Table 
1  both  adherends  were  glass.  Preparation  of  the  glass  cover  slip  surfaces 
is  less  well  standardised.  The  fundamental  problem  is  removal  of  oxides 
of  low  surface  energy.  Chemical  attack  by  common  oxidising  agents,  e.g. 
nitric  acid  and  hydrogen  peroxide  certainly  remove  a  large  fraction  of 
these  oxides.  So  too  does  exposure  to  the  oxidising,  region  of  a  naked 
flame.  Sand  blasting,  in  order  to  roughen  the  surface  and  promote 
mechanical  linkages,  and  application  of  coupling  agents  (silanes)  have 
also  been  recommended.  However,  no  combination  of  these  treatments  has 
yet  produced  a  joint  which,  when  it  eventually  debonds  after  prolonged 
exposure  to  water  in  the  absence  of  externally  applied  stress,  leaves 
fracture  surfaces  which  do  not  include  large  areas  of  glass/adhesive 
interface . 

In  joints  manufactured  using  liquid  resin,  a  reproducible  glue  line 
thickness  was  achieved  by  simply  allowing  the  cover  slip  to  sink  under  its 
own  weight.  In  joints  made  using  solid  adhesive  film,  the  glue  line 
thickness  was  controlled  by  clamping  during  the  cure  process. 

Adjacent  rings  of  the  same  colour  (black  or  white)  ate  loci  of  points  for 
which  the  optical  path  length,  in  the  space  between  cover  slip  and  optical 
flat,  differs  by  one  wavelength.  By  the  same  token,  a  displacement  in  the 
pattern  of  Newton's  rings  by  an  amount  equal  to  one  ring  width  corresponds 
to  a  change  in  path  length  equal  to  one  wavelength.  By  observing  changes 
In  the  number  of  rings  between  fixed  markers,  such  as  entrapped  air 
bubbles,  displacements  normal  to  the  joint  during  water  uptake  can  be 
measured  to  an  accuracy  of  A/4.  If  required,  displacements  which  are  at 
least  as  small  as  A/10  can  be  measured  by  super  impos  i  t  ion  of  images  in 
order  to  create  Moire  patterns. 

The  water  concentration,  and  hence  the  swelling  associated  with  water 
uptake,  eventually  saturates  and  thereafter  a  shoulder  separating  fully 
saturated  from  less  than  fully  saturated  resin,  progressively  moves 
inwards  from  the  rim  of  the  joint.  Tills  is  shown  schematically  in  Figure 
5c.  Figure  6  shows  a  sequence  of  images  obtained,  showing  the  development 
of  swelling  in  a  joint  manufactured  between  anodized  titanium  and  a  glass 
cover  slip  using  adhesive  B  and  immersed  in  distilled  water  at 
approximately  80°C. 


To  maintain  contact  with  the  outer  ring  of  uniformly  swollen  adhesive 
(Figure  5c),  the  adherends  would  each  need  to  bend  with  curvature  that  is 
opposite  to  that  inside  the  shoulder.  Failure  to  adopt  such  "S"  wise 
bending  manifests  itself  as  an  interfacial  crack,  the  occurrence  of  which 
is  revealed  as  patches  of  bright  contrast  at  the  rim  in  Figure  6. 


The  positions  of  well-defined  points  in  the  patterns  of  Newton's  rings  may 
be  measured  as  functions  of  time.  Figure  7  shows  a  plot  of  the  shoulder 
defining  the  extent  of  water  saturation,  and  of  the  crack  edge  as 
functions  of  (time)  for  an  aluminium  joint  manufactured  with  adhesive  A 
and  immersed  in  distilled  water  at  20  C. 


Figure  8  shows  a  plot  if  the  shoulder  defining  the  extent  of,  water  ingress 
form  the  experiment  reported  in  Figure  6. 

Figure  9  shows  a  similar  plot  for  a  specimen  manufactured  using  adhesive  B 
but  with  a  stainless  steel  substrate  and  for  an  immersion  water 
temperature  of  80°C.  From  the  marked  similarity  in  the  results  for 
anodized  titanium  and  stainless  steel  it  is  concluded  that  the  deformation 
is  wholly  within  the  adhesive  and  is  independent  of  the  adherend.  Moire 
patterns,  generated  by  superimposition  of  interference  patterns  before  and 
after  swelling  were  used  to  obtain  the  data  in  Figure  10.  These  specimens 
were  manufactured  from  adhesives  C  and  D  and  measurements  were  made  at  a 
range  of  temperatures  between  ambient  and  100°C. 

If  migration  of  water  occurs  at  low  water  concentrations  and  is  governed 
by  Fick's  law,  plots  of  migration  distance  as  a  function  of  (time)^" 
should  approximate  to  straight  lines.  The  data  points  for  specimens 
manufactured  with  adhesives  A  and  B  certainly  do  not  approximate  to 
straight  lines.  However  those  manufactured  with  adhesives  C  and  D  and 
which  were  studied  using  the  more  precise  Moire  measuring  technique  do 
suggest  Fickian  behaviour  although  the  slopes  of  the  graphs  indicate 
activation  energies  that  are  much  higher  than  normally  expected  for  epoxy 
resins. 

To  draw  a  comparison  of  this  apparent  Fickian  behaviour  in  the  specimens 
with  the  theory,  we  can  solve  Fick's  law  in  cylindrical  polar  coordinates, 
assuming  axial  symmetry: 

<)£  =  t2C  +  ±  bC 

dT  =  iR  R  3R 

where  T  =  Dt/a  ,  R  =  r/a,  C  =  c/cgAT  with  D  a  constant  diffusion 
coefficient. 
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By  using  a  suitable  substitution  this  equation  is  Bessel's  of  order  zero, 
and  a  solution  may  be  found.  For  small  values  of  T,  the  infinite  series  is 
slow  to  converge  and  computation  is  impractical.  As  T  increases  however, 
it  is  possible  to  obtain  series  which  converge  quite  rapidly.  Such  a 
series  has  been  used  to  obtain  the  theoretical  curve  of  Figure  11  for  a 
Redux  312/5  epoxy  at  62°C  where  the  time  is  334  hours.  It  was  found 
experimentally  that  full  saturation  was  represented  by  14  fringes  and  the 
experimental  results  have  been  included  in  the  Figure  11. 
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3. 


SELF  STRESSED  FRACTURE  MECHANICS  TEST-1MECE 


3.1  THE  EEASTIC  DISC  IN  PLANE  STRAIN 

The  disc  consists  of  a  circle  of  aluminium  which  in  plane  polar  co¬ 
ordinates  is  defined  by  the  relations  f)<r<a  and  0<P<2n.  It  is  stressed  bv 
removing  a  sector  fKr<a,  2m-fl<0< 2 it  and  making  an  adhesive  joint  wlien  the 
pap  has  been  closed.  Assuming  that  plane  stress  conditions  are  met,  it  is 
necessary  to  calculate  the  stress  intensity  factor,  Kj ,  and  the  formation 
energy,  W,  for  the  crack  c<r<a,  0  =  0.  Hi  is  has  been  done  by 
superimposing  the  solutions  of  the  resulting  stress  with  those  of  K.  and  U 
for  a  traction-free  circumference  and  a  symmetrical  pressure  loading. 


The  Airy  function  given  by  Mitchell  (1899),  Eshclby  (1966) 

'F(r,0)  =  aQ  +  bQ  lnr  +  cQr2  +  dQr2lnr  3.1 

may  be  used  to  find  the  stress  and  displacement  fields: 

orr(r,0)  =  bor-2  +  2cq  +  dQ(l  +  2  lnr)  3.2 

arg(r,0)  =0  3.3 

o00(r,0)  =  -bQr-2  +  2cq  +  dQ(3  +  2  lnr)  3.4 

ur(r,Q)  =  E"1 {-b0(l+v)r~1+2co(l-v)r+2d0(l-v)r  lnr-d0(l+v)r}  3.5 


and 


u  q ( r , 0 )  =  4E-1  dQrO  3.6 

Since  the  solution  is  required  to  satisfy  the  conditions 

(i)  orr(a,0)  =  orO(a,0)  =  0,  O«0>;2n 

(ii)  Uq(c,  =  r!2  an^ 

(iii)  u  and  iiq  are  bounded  at  r  =  0,  we  see  that  bQ  =  0,dQ  =  E°/8n 

and  cQ  =  ~(EC/8n)  ln(a/e). 
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Therefore 


T  (r,0  ) 


and 


TO  r*"  ln(r/a  /e) 
8 


oD0(r,O)  -  Bp  { 1  +  I n ( r / a ) } 

An 


3.  7 


3.8 


It  follows  that  we  may  wrfte  o(r)  =  pQf(r/a) 
where  pQ  =  K)/4n 
and  f(p)  =  1  +  Inp 

To  find  Kj  and  W  the  dimensionless  quantities  p  =  r/n  (already  introduced) 
may  be  used 

spp(p,0)  =  orr(r,0),  spQ  (p ,  0 )  =  orf)(r,0),  s0Q(p,O)  =  op0(r,n)  3.9 

P  o  Po  Pq 

u(p,0)  =  F.ur(r,0)  and  v(p,0)  =  F.u^Cr.O) 

Pq  O+v  )a  pQ  ( l+v  )a 

then,  by  making  use  of  a  symmetry  argument,  it  is  easily  shown  that  the 
problem  to  find  and  W  may  be  stated  as  follows: 

Solve  the  dimensionless,  plane  stress  equations  of  elasticity  in  the  semi¬ 
disc  0<p<l,  0<0<it  subject  to  the  conditions 

i)  The  stresses  and  displacements  are  bounded  at  r  =  0 

ii)  spp(l,0)  =  spO(l,0)  =  0,  0<0<n 

iii )  sp0 (P»0)  =  s  0(p,n)  =  0,o<p<l 

iv)  v(p,rt)  =  0,  0<p<l 

v)  v(p,0)  =  0,  0<p<c 

vi)  s00(p,O)  =  -  (1  +  log  p),  c<p<l 

vii)  limit  f1  v  ( p ,  0 )  <“> 

p  +  i  dp 

By  slightly  modifying  the  results  in  section  2  of  Tweed  and  Rooko  (1973) 
it  can  bo  shown  that  there  is  a  solution  of  the  equations  of  elasticity, 
in  the  given  region,  which  satisfies  the  conditions  (i)  through  (v)  and  is 
such  that 
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see  (p,0>  =  1 


P(t)  M(p,t)dt, 
/[(l-t)(t-c)J 


0<p<  1 


3.  10 


and 


v(p , 0)  =  2 


P(t)  dt 


I+v  J  /((l-t)(t-c)] 
c 


c<p<  1 


3.11 


where  P(t)  is  an  arbitrary  function  of  t 
and 

M(p,t)  =  1  +  (t2  -  l)2  +  t(t2  -  1)  -  _t_  +  t  -  _1_ 

t-p  t(l  -  pt)3  (1  -  pt)2  (l  -  pt)  t 


3.12 


It  follows  that  conditions  (vi)  and  (vii)  will  be  satisfied  also  if  P(t) 
is  a  solution  of  the  singular  integral  equation. 

1 

I  fp(OM(p.t)dt  =  -  (1  +  log  p),  .13 

n  J  /[ ( 1 - t ) ( t -c ) ] 

c 


with  subsidiary  condition 

P  ( 1 )  =  0 


3.  14 


The  stress  intensity  factor  Kj  and  the  crack  formation  energy  W  are 
defined  by  the  equations 

Kj  =  limit  /2it(r-ac)  F,  3up(r,0)  3.15 

r-*-ac+  2  3r 

and 

W  =  J  p0 [ 1  +  log (r /a) ]  ug(r,0)dr  3.16 

ac 

respectively.  Therefore,  by  3.10  and  3.11  we  see  that 

KI  =  Po  ^  2)1/1  P(c)  3-17 

i  -  c 
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and 


1 


W  =  -  2po^a^  j  tP(t)logt  dt 
✓  l(l-t)(t-c)] 


3.  18 


Let  K0  and  W  be  the  stress  Intensity  factor  and  fornation  energy 
respectively  of  a  crack  of  length  2 (a-ac)  in  an  infinite  elastic  solid, 
when  the  crack  is  opened  by  a  uniform  pressure  p  ,  then  after  Sneddon  and 
Lowengrub  (1969) 


=  P_  /  Tia(l-c) 


3. 19 


and 


n  -  2  2  .2 

W0  =  up0  a  (1-c) 


3.20 


Therefore 


Kj.  =  /2_  P(c) 
K  1-c 


3.21 


and 


,1 

t  P(t)  log  t  dt 


W  =  -2 

Wc  7C(1-C)  2  J/[(l-t)(t-c)] 

c 


22 


Also  of  interest  is  the  strain  energy  release  rat which  is 
defined  by 


fjU)  -  OW  -  1  Kj' 


di  E 

where  £  =  a(l-c)  is  the  crack  length. 
If  we  define 

1  Ko2  =  Po2 

r  12 


tlicn 


=  h. 


Jo  Kq‘  (1-c) 


2 _  P*(c) 

2 


3,2  3 


3.2 


Numerical  solutions  to  equation  3.5  are  found  by  the  Caviss-Chobyshev 
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quadrature  technique,  Erdogan  and  Gupta  (1972),  and  the  ratio  Kj/Ko  using 
Chebyshev-Lagrange  interpolation.  The  Gaussian  quadrature  formula  is  used 
to  determine  W/W  . 


3.2  EXPERIMENTATION 

A  pilot  scheme  was  executed  to  investigate  the  parameters  controlling  the 
design  of  the  self  stressed  fracture  mechanics  test-piece.  SIC  aluminium 
discs,  3"  (7.6  cm)  diameter,  1/16"  (0.16  cm)  thick  with  a  5°  sector 
removed,  provided  the  first  test-pieces.  In  practice  the  sector  in  the 
pilot  test-pieces  was  closed  by  deforming  the  disc  into  a  cone  which, 
after  joining  the  cut  faces,  was  flattened  by  clamping  between  a  pair  of 
plattens.  The  theory  demands  boundary  conditions  that  naturally  set  the 
shear  stresses  and  radial  component  of  stress  at  the  free  surface  to  zero, 
and  that  the  sector  is  closed  up  simply  by  the  90  component  of  stress. 

These  original  test-pieces  were  found  to  be  stressed  beyond  the  elastic 
limit  in  some  areas  of  the  discs  and  further  aggravation  came  from 
buckling.  Mitchell  and  Head  (1961)  calculated  that  a  disci inated  plate 
will  buckle  if  the  angle  -ft  buckle  exceec*s  ^  (thickncss/radius)".  For  the 
original  disc,  =  1°20',  therefore  buckling  presents  a  problem.  It 

became  necessary  to  optimise  the  three  parameters: 

(i)  sector  cut  out  angle 

(ii)  disc  radius 

(iii)  disc  thickness. 

Variations  in  (i)  and  (ii)  lead  to  difficulties  in  preparing  a  sound  joint 
and  thus,  the  thickness  was  increased  to  3/16"  (0.48  cm)  which  removes  all 
chance  of  buckling.  The  discs  could  not  now  be  manufactured  by  punching, 
each  having  to  be  machined,  and  although  this  will  load  to  a  more 
expensive  specimen,  it  should  bo  borne  in  mind  that  the  test-piece  can  be 
used  many  times  as  it  functions  wholly  within  the  elastic  regime. 

Closing  the  gap  to  allow  a  joint  to  be  made  at  the  cut  has  been  very 
carefully  examined-  A  custom  built  chain  wrench  closes  the  gap  too 
quickly,  and  a  range  of  differing  gaps  were  tried  to  solve  the  problem.  To 
compensate  for  the  pinching  effect  a  parallel  sided  cut  was  effected  using 
a  range  of  gap  widths.  This  reversed  the  pinching  effect  in  that  closure 
occurred  first  at  the  circumference.  A  series  of  tests  was  also  carried 
out  where  the  elastic  limit  was  exceeded  and  the  amount  of  plastic 
deformation  was  monitored.  The  data  from  these  results  allowed  the 
dimensions  of  a  wedge  to  be  determined  which  would  maximize  the  amount  of 
elastic  strain  energy  stored  in  the  disc,  and  allow  the  two  faces  to  close 
up  evenly  within  the  no  buckling  requirement.  The  wedge  dimensions  were 
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set  at  1.1  mm  (0.046")  at  the  circumference  and  0.75  mm  (0.032")  at  the 
Inner  edge,  giving  a  sector  angle  approximately  equal  to  1°. 

The  cleaning  process  carried  ojt  immediately  prior  to  the  adhesive  joint 
formation  was  as  follows: 


TIME  TEMPERATURE 

(mins)  (°C) 


1. 

Chlorinated  solvent 
(carbon  tetrachloride) 

5 

ambient 

2. 

Non-si l icated  alkaline  cleaner 
(0.1  molar  sodium  hydroxide) 

10-15 

6  3 

3. 

Tap  water  immersion 

O 

ambient 

4. 

Tap  water  immersion 

5 

38-60 

5. 

Deoxidiser  (170  ml  nitric  acid 

30  ml  hydrofluoric  acid,  800  ml  water) 

10-15 

amb i ent 

6. 

Tap  water  immersion 

a 

ambient 

7. 

Deionised  water  spray 

- 

ambient 

8. 

Phosphoric  acid  9-11%  by  wt  in 
deionised  water 

- 

ambient 

Apply  voltage  for  1  minute.  Step  to  10  volts 
after  2-5  minutes.  Maintain  for  20-25  minutes 

9. 

Tap  water  spray 

10-15 

ambient 

10. 

Deionised  water  spray 

2 

amb i on  t 

11. 

Dry  in  air 

30 

80 

12. 

Cool  in  air 

120 

ambient 

13. 

Prime  within  120  min 

- 

- 

14. 

Dry  in  air 

30 

ambient 

15. 

Dry  in  air 

60 

123 
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The  cleaned  aluminium  test  discs  had  the  epoxy  resin  strip  or  layer  of 
paste  positioned  between  the  opposing  faces  of  the  wedge  cut  out,  with 
restraining  plates  on  each  side  of  the  discs  to  maintain  flatness  during 
the  closing  operation.  Standardisation  was  achieved  using  a  torque  wrench 
set  at  16  ft-lb.  Immediately  on  closure,  the  restraining  plates  were 
removed  and  the  prepared  joints  gelled  and  cured.  The  manufacturer's  data 
sheets  were  strictly  followed  for  both  cure  times  and  temperatures. 
Examination  of  the  joints  revealed  reasonable  reproduction  of  the  bond 
line  thickness,  approximately  the  thickness  of  the  resin  strips  used, 
0.25  mm.  The  proprietary  adhesive  resin  film  most  extensively  employed  was 
FM73M,  but  others  which  were  tested  in  the  experimental  programme  were 
Redux312/5,  and  DLS468/DLS476. 

The  length  of  crack  when  present,  was  measured  and  the  disc  was  then 
immersed  in  boiling  water  for  varying  periods  of  tine.  After  each 
immersion  the  length  of  crack  was  measured.  This  process  was  repeated 
until  the  crack  had  propagated  the  full  length  of  the  glue  line,  with 
frequent  water  changing  to  avoid  contamination.  On  examining  the  crack  it 
was  always  found  that  the  crack  length  was  different  when  viewed  from  the 
two  sides  of  the  disc.  For  the  purpose  of  evaluation,  the  average  value  of 
the  crack  length  read  from  both  sides  was  noted.  Also,  near  the  crack  tip 
it  became  increasingly  difficult  to  distinguish  where  the  crack  actually 
finished,  a  result  of  the  hairline  nature  of  the  crack.  When  the  crack 
length  was  apparently  shorter  in  the  glue  line  after  further  boiling,  the 
length  of  crack  was  taken  to  be  the  distance  to  the  end  of  the  last  crack. 
With  some  discs  that  did  not  initially  crack  on  release  of  the  clamp  an 
attempt  was  made  to  induce  a  starter  crack  in  the  glue  line  by  cither 
sawing  a  small  nick  at  the  periphery  of  the  joint,  or  1  •  immersing  the 
discs  in  liquid  nitrogen  to  increase  the  stress  field  in  the  disc  by 
contraction  of  the  aluminium.  Unfortunately,  neither  method  proved  to  be 
successful. 

A  series  of  experiments  were  performed  using  a  perspex  disc  of  similar 
dimensions  to  the  aluminium  specimens.  The  object  here  was  to  use  a 
photoelastic  method  to  study  the  stress  field  during  crack  growth. 
Initiation  of  the  crack  has  proved  difficult  and  satisfactory  monitoring 
of  the  stress  field  during  propagation  is  presently  under  study. 

Two  undergraduate  students,  Cilbert  and  Longden,  working  hero  in  Bristol 
have  become  particularly  Interested  in  the  experimental  evaluation  of  the 
fracture  mechanics  test-piece  and  are  taking  this  piece  of  interesting 
research  a  stage  further. 


27 


3.3 


RESULTS  AND  DISCUSSION 


On  release  of  the  chain-wrench  damp  from  around  the  stressed  disc,  crack 
lengths  vary  according  to  the  epoxy  resin  adhesive  formulation  used  to 
make  the  joint.  The  crack  always  propagates  the  full  glue- line  length 
(1/a  =  1)  Mr  resin  Redux  312/5,  whereas  for  F1173M  no  crack  is  present  in 
the  joint  tl/a  =  0).  The  1 /a  values  obtained  for  IH.S46S /DLS4  7  6  are  0,  0, 
0,  0.40  and  0.82. 

Correlation  of  crack  length  with  independently  measured  values  of  tjjg  is 

required  in  order  to  generate  calibration  curves  for  a  given  geometry  of 

sel  f -s  t  ressed  disc.  Available  data  for  the  dry  adhesives  show 

variations,  e.g.  FM73M  has  values  of  2.38  k.J/M-  and  3.67  kd/M*"  for  two 

different  batches.  Redux  312/5  and  DLS468/DI  S476  both  have  Cj  ,  values  ot 
a  a  ■'I1, 

2.8  kJ/M"-.  Using  as  calibration  points,  the  tj  j  values  for  dry  FM73M  and 

Redux  312/5  and  the  stable  crack  lengths  measured  after  mechanically 

introducing  started  cracks,  two  values  of  1 /a  measured  for  DI.S468/DLS 47b 

indicate  that  (  tjip)max  is  3.65  +  0.55  kJ/M“. 

Conflict  evidently  exists  between  the  published  data  and  the  results 
obtained  here;  with  DLS468/DLS476  our  self-stressed  discs  were  lound  to 
crack  partially  or  not  at  all.  With  Redux  312/5,  the  bonds  crack  over  the 
full  length  of  the  glue  line.  With  FM73M  no  cracks  are  found  suggest i up. 
that  is  greater  than  2.8  kJ.M“.  With  a  view  to  clearing,  up  this 

dilemma,  independent  measurements  of  on  the  particular  batches  of 

adhesive  films  used  here  are  being  made  at  the  University  of  Cambridge. 

Results  presenting  the  variation  of  Cjj^,  expressed  as  a  fraction  of 
(  yi(_-)n.lx  with  water  uptake  by  FM73M  are  shown  in  Figures  12  to  1  The 
specimens  were  immersed  in  boiling  distilled  water.  The  tine  for  the 
starter  crack  to  appear  varied  from  5  to  130  mins.  Examination  of  the 
joints  after  total  failure  show  that  propagation  of  the  cracks  occur 
through  the  resin  and  definitely  not  at  the  resin/aluminium  interface. 
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4. 


PHOTOELASTIC  IMAGING  OF  DEFECTS 


Andrews  and  Wallis  (1977)  were  first  to  demonstrate  that  a  light 
emitting  diode(LF.D)  can  be  used  to  illuminate  the  stress  waves  travelling 
in  a  visualizing  block  when  viewed  through  crossed  polars.  The  method  was 
applied  to  the  detection  of  defects  in  opaque  solids  by  Hardy,  Turner  and 
Asiibee  (1978).  There  is  a  strong  reflected  component,  mode  converted, 
from  a  defect  which  can  be  observed  some  distance  from  the  flaw,  after  the 
stress  wave  has  travelled  from  the  specimen  into  the  transparent  block. 
T!i  i  s  wave  interacts  with  direct  waves  from  the  transducer  to  give 
interference  patterns  for  which  a  computer  model  has  been  devised  by 
Kitson,  Low  and  Turner.  All  the  longitudinal  and  shear  waves  emanating 
from  the  transducer,  the  interface  and  the  defect,  wit li  and  without  mode 
conversion,  are  investigated. 


4.1 


COMPUTER  MODEL 


Figure  16  shows  the  representation  of  the  model,  witli  a  series  of  radiating 
point  sources  depicting  the  transducer.  An  array  of  points  P(x,v)  forms 
the  visualized  zone  and  was  originally  set  at  31  x  31  to  conform  with 
computer  graphics.  It  is  assumed  that  the  defect  is  perfectly  round,  that 
longitudinal  (L)  and  shear  (S)  waves  emanate  from  the  transducer  point 
sources,  and  that  on  reflection  and  refraction  both  mode  converted  and 
unconverted  waves  emerge.  The  following  relationships  hold: 


sin  <x  , 


V 


S 


s  i  n  V  =  s  i  n  <5 


where  and  V,,  are  the  velocities  in  the  specimen  and  11,.  the  shear 
velocity  in  the  visualizing  quartz.  Note  that  the  subscripts  arc  inter 
changeable  for  the  incident  shear  wave  reflected  as  longitudinal  and  for  a 
refracted  longitudinal  wave.  Relative  intensities  may  be  evaluated  using 
the  Zoeppritz  equations  and  an  attenuation  factor  of  l//x  is  included  for 
tlie  wave  amplitude  as  it  travels  away  from  the  source.  Hence  we  may  find 
the  exact  path  which  any  particular  ray  travels  together  with  its 
amplitude  at  a  point.  Each  flash  of  the  LED  is  represented  by  a  sine  wave, 
singly  for  single  flash  and  several  cycles  for  multi-flash. 


The  complete  picture  is  built  up  as  follows:  each  typo  of  wave 
[longitudinal/shear  (LS),  long/long/shear  (LI,S),  etc.]  is  considered  in 
turn  for  each  transducer  point  source.  The  program  calculates  (1)  the 
patli  between  the  transducer  point  and  a  point  in  the  rectangular 
visualized  array  (subject  to  the  path  being  possible,  (ii)  when  the  sound 
pulse  arrives  at  the  visualizing  point,  and  (iii)  the  amplitude.  The 
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amplitude  is  then  multiplied  by  the  value  of  the  illumination  waveform  at 
that  specific  time,  and  repeated  for  all  the  points  in  the  rectangle. 
Results  for  all  the  transducer  elements  are  superimposed,  and  finally  the 
whole  process  repeated  for  each  wave  type  in  order  the  simulate  the 
complete  interference  pattern.  The  final  array  is  normalized  to  give  the 
negative  or  compressive  regions  the  lighter  shades  of  the  computer 
simulation  and  the  positive  or  tensile  regions  the  darker  areas:  a 
normalization  analogous  to  the  quarter  wave  plate. 

Two  main  calculation  procedures  are  used  in  the  program,  RF.ITATH  and 
DIRCALC.  REIT  ATI),  the  flow  diagram  is  shown  in  Figure  17,  deals  with  the 
patterns  resulting  from  the  wavetypes  which  are  reflected  from  the  defect, 
whilst  DIRCALC  handles  those  direct  from  the  transducer.  Any  wavetypes 
can  be  included,  but  it  is  found  experimentally  that  many  are  of  little 
importance  and  only  LLS,  LS,  and  SSS  need  to  be  considered.  It  should  be 
remembered  here  that  the  experimental  arrangement  of  the  photoelastic 
method  helps  enormously  to  reduce  the  number  of  wavetypes  which  contribute 
to  the  visualization.  In  the  direction  of  principal  stress,  i.e.  the 
direction  of  the  stress  wave  propagation,  45  degrees  to  the  polar  axes, 
only  longitudinal  waves  are  visualized,  and  at  45  degrees  to  this 
direction,  along  the  polar  axes,  only  transverse  waves  can  he  seen.  It  is 
in  this  direction,  along  the  polar  axes  that  we  are  primarily  concerned  as 
it  is  here  that  most  information  is  displayed. 

The  first  attempts  to  display  the  visualized  array  by  computer  graphics 
took  advantage  of  SYMAP  (SYMAP  is  a  computer  program  for  producing  maps 
which  graphically  depict  spacially  disposed  quantitative  and  qualitative 
information.  SYMAP  was  designed  and  developed  by  the  Laboratory  for 
Computer  Graphics  at  Harvard  University,  USA]  and  then  SYMVU  was  employed 
to  give  a  topographical  representation  of  the  data.  The  former  program 
has  the  advantage  that  direct  comparisons  could  be  made  of  the 
photographic  results,  but  suffered  from  having  only  ten  levels  of  shading, 
whereas  the  latter  shows  the  dislocations  very  well.  Figure  IS  shows  the 
contour  map  where  the  maximum  pressure  is  represented  by  a  height  of  3 
ins,  with  a  minimum  of  0  ins.  The  viewing  point  is  at  an  altitude  of  30 
ins  with  an  azimuth  of  335  degrees.  Dislocations  in  this  context  refer  to 
the  line  of  discontinuities  in  wavefronts  which  line  runs  in  the  direction 
of  the  stress  propagation.  A  new  program  was  written  therefore,  following 
the  principles  of  SY'MAP,  but  with  facility  for  100  shading,:'..  Some  time 
was  spent  experimenting  with  different  overprinted  characters  in  order  to 
get  a  set  of  evenly  varying  shades.  It  was  found  that  a  sub-routine 
designed  to  interpolate  the  brightness  between  the  array  points  gave  a 
smoother  picture  when  the  visualization  density  map  was  plotted  on  the 
Line  printer. 
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4.2 


VISUALIZED  SOUND  EXPERIMENTS 


Figure  19  shows  a  schematic  diagram  of  the  apparatus.  The  LED  sends  a 
single  pulse  (single  flash)  or  flashes  a  number  of  times  in 
synchronisation  with  the  ultrasonic  frequency  (multi-flash)  to  illuminate 
the  stress  field  in  the  quartz  visualizing  block.  The  birefringence 
exhibited  by  quartz  when  it  is  stressed,  results  in  a  snail  component  of 
light  being  propagated  perpendicular  to  the  axis  of  the  polarizer  to  give 
elliptically  polarised  light.  This  contrasts  with  transmission  of  the 
illuminating  pulses  through  uns t ressed ^regions  of  the  quartz  since  here 
the  light  remains  plane  polarised.  The  analyser,  set  at  90  degrees  to  the 
polarizer,  permits  transmission  of  the  small  component  only  so  that  the 
stressed  regions  are  revealed  as  bright  fringes.  Bright  fringes  arise  in 
both  tensile  and  compressive  regions,  and,  by  suitably  setting  the  quarter 
wave  plate,  all  the  light  from  the  most  highly  compressed  region  can  be 
eliminated.  Thus  some  light  can  be  seen  for  unstressed  regions  with  the 
tensile  areas  brighter  still,  the  whole  picture  being  built  up 
stroboscopically  at  repetition  rates  up  to  3  kHz.  Full  details  of  the 
experimental  approach  are  given  in  Hardy,  Turner  and  Ashbee  (1978). 

The  experiments  which  have  been  carried  out  in  parallel  with  the  computer 
modelling  have  been  mostly  single  flash,  i.e.  a  single  LED  pulse  for 
illumination.  Single  flash  images  are  difficult  to  record  because  they 
are  so  dim,  but  it  is  necessary  to  identify  the  different  wavetypes.  A 
series  of  experiments  were  conducted  with  no  specimen  in  place  in  order  to 
investigate  the  transducer  for  accurate  computer  modelling.  For  these, 
three  flashes  were  used.  As  the  shear  waves  travel  more  slowly  than  the 
longitudinal  waves,  the  shear  fringes  arc  closer  together:  a  fact  which 
with  their  position  simplifies  the  identification  of  any  particular  wave. 
To  investigate  the  transducer  beam,  a  series  of  photographs  was  taken 
varying  the  polar  axes,  thus  compensating  for  the  angular  sensitivity  of 
the  optical  arrangement.  The  angular  spread  of  the  beam  was  found  to  be 
approximately  90  degrees. 

A  number  of  different  transducers  have  been  tried  including  two  commercial 
designs,  but  the  results  have  been  obtained  witli  a  transducer  designed  and 
built  in  the  laboratory.  It  consists  essentially  of  a  brass  plate  0.9 
ins.  (2.2  cm)  diameter  and  5^/4  thick  on  to  which  a  0.4  in.  (1.0  cm) 
diameter  piezo-elcctric  disc  is  soldered.  However,  the  imago  is  still 
complicated  by  secondary  waves  arising  from  reflections  within  the  brass 
and  di sc. 

After  the  initial  experiments  it  appeared  that  two  wavetypes  were 
contributing  to  interference  patterns  seen  in  the  visualized  region. 
Firstly,  the  wave  reflected  from  the  defect.  This  is  mode  converted  blit 
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crosses  the  interface  between  specimen  and  visualizing  block  without  a 
further  conversion  (LSS).  Secondly,  the  direct  wave  originating  as  shear 
at  the  transducer  and  again  crossing  the  interface  without  mode  conversion 
(SS).  The  experiments  show,  however,  that  the  prominent  shear  waves  leave 
the  spec  i  men /qua  r  t  z  interface  at  roughly  the  same  time  as  the  primary 
longitudin  ils.  To  investigate  this  further,  the  positions  of  the  primary 
longitudinal  and  shear  waves  were  plotted  for  a  series  of  different  delay 
times.  Different  specimen  thicknesses  were  used  for  the  measurements  as 
well  as  readings  taken  without  a  specimen.  The  results  shown  in  Figure  20 
strongly  suggest  that  the  shear  waves  arise  by  way  of  mode  conversion  at 
the  interface,  having  crossed  the  specimen  as  longitudinal  waves.  Shear 
waves  generated  from  the  transducer  are  apparently  unimportant.  Different 
effects  of  the  quarter  wave  plate  on  the  two  wavetypes  account  for  the 
incorrect  positioning  of  the  specimen/quartz,  interface  in  that  Figure. 

To  establish  a  baseline  from  which  to  judge  images  from  resin  specimens  a 
range  of  steel  specimens  have  been  produced  with  thicknesses  up  to  1  1/4 
ins  (31.2  mm)  each  containing  a  transverse  hole  'defect'  having  diameters 
in  the  range  of  0.04  ins  (1  mm)  to  0.4  ins  (10  mm).  A  35  mm.  Nikon  F2 
photomic  camera  is  used  to  record  the  results  with  an  eight  minute 
exposure  on  ASA  400  film  for  single  flash  at  different  delay  tiros. 

The  delay  time  is  measured  from  the  original  ultrasonic  pulse  to  the  LED 
flash,  and  is  the  time  taken  for  the  sound  to  travel  from,  the  transducer 
to  the  visualized  position.  The  TTL  control  enables  this  to  he  accurately 
measured,  since  the  delay  monostable  and  thyristor  fired  ultrasonic  wave 
are  triggered  by  the  same  pulse.  A  digital  readout  of  the  delay  is  used 
which  can  be  checked  against  the  tine  base  of  the  oscilloscope. 


4.3  RESULTS  AND  DISCUSSION 

The  results  are  presented  in  a  graphical  form  which  show  both  experimental 
position  and  the  computer  simulation.  Figures  21  and  22  show  the  results 
for  two  delays,  8.8  and  13.4  s  respectively  for  steel  specimens.  The 
dotted  lines  give  the  position  of  the  computer  prediction  and  the  hard 
lines  have  been  traced  from  the  photograph.  They  show  that  the  program 
correctly  predicts  the  wavefront  spread  for  both  the  LSS  and  F.S  waves.  It 
is  the  former  which  represents  'information'  from  the  defect.  Differences 
in  the  lengths  of  the  wavefronts  can  be  attributed  to  the  differences  in 
the  transmission  coef fecients.  The  latter  are  complicated  by  the  acoustic 
coupling  at  the  spec i men /q uar t z  interface  by  a  thin  film  of  water. 
Broadening  of  the  lower  LS  wave  is  caused  by  the  superposition  of  the  LL5 
wave  i.e.  the  reflection  from  the  hole  without  mode  conversion. 


The  experimental  delay  times  against  the  theoretical  time  required  to 
align  the  LS  waves  is  shown  in  Figure  23.  The  gradient  of  unity  within  3Z 
demonstrates  further  that  the  program  correctly  predicts  the  LS  position. 

A  sequence  of  tests  was  undertaken  with  thicker  specimens  1  1/4  ins  (31.2 
mm)  containing  a  hole  0.3  ins  (8  mm)diameter,  where  the  distance  between 
the  hole  centre  and  the  transducer,  A,  was  varied.  Figures  24,  25  and  26 
give  the  results  for  A  =  0.27  ins  (7  mm),  0.61  ins  (15.6  mm)  and  0.95  ins 
(24.2  mm)  respectively. 

The  above  results  refer  to  single  flash  which  was  regarded  as  the  best 
approach  for  wave  indentif ication,  whereas  the  original  approach  began 
with  multiflash.  More  information  is  carried  in  the  multiflash  pictures 
as  the  close  up  of  the  interference  patterns  from  0.04  in  (1  min)  radius 
'defect'  shows  in  Figure  27.  The  multiflnsh  simulation  of  Figure  28 
illustrates  the  interaction  of  the  LS  wave  and  the  LSS  wave  reflected  from 
a  hole  of  0.04  ins  (1  mm)  radius. 

In  composite  materials  it  has  been  possible  to  calculate  both  longitudinal 
and  shear  wave  velocities  by  measuring  the  delays  across  a  known  specimen 
thickness.  Velocity  variations  have  been  detected  during  water  uptake  and 
are  attributed  to  resin  plastication  and  to  post  curing.  Experiments  on 
composite  laminates  containing  holes  have  been  less  successful. 
Attenuation  of  the  wavetypes  is  more  marked,  and  the  shear  wave  rapidly 
loses  energy  as  it  proceeds  across  a  series  of  interfaces.  The  problem  is 
aggravated  by  higher  frequencies,  resolution  of  the  fringes  being 
impossible  at  frequencies  in  excess  of  ~ 10  MHz. 
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General  view  of  ambient  temperature  apparatus. 


Hi ph  temperature  (<100°C)  apparatus. 
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epoxy 


STAINLESS  STEEL  FM  1000  ~80  C 


Figure  9.  Diffusion  rate  for  an  FM  1000  adhesive  joint  between  stainless  steel  and  glass  cover  slip. 
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Computer  model  representation. 


Figure  19.  Schematic  diagram  of  photoelastic  system. 


tit  A 


3  x  3  cm 


*<i 


(x,y)  =  (1-50, 2-00)  cm 


Photograph  - Computer  prediction 

Photograph- Printout  Comparison 
Delay  =  8*  Bps 


Figure  21. 
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Figure  22. 


Computer  simulation. 
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Figure  24.  Computer  simulation. 
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Figure  25.  Computer  simulation. 
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Figure  26. 


Computer  simulation. 


